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Abstract. We use the clathrate hydrate trapping theory to calculate the enrichments in O, N, S, Xe, Ar and 
Kr compared to solar in Saturn's atmosphere. For this, we calibrate our calculations using two different carbon 
abundance determinations that cover the domain of measurements published in the last decades: one derived 
from the NASA Kuiper Airborne Observatory measurements and the other obtained from the Cassini spacecraft 
observations. We show that these two different carbon abundances imply quite a different minimum heavy element 
content for Saturn. Using the Kuiper Airborne Observatory measurement for calibration, the amount of ices 
accreted by Saturn is found to be consistent with current interior models of this planet. On the other hand, using 
the Cassini measurement for calibration leads to an ice content in the planet's envelope which is higher than 
the one derived from the interior models. In this latter case, reconciling the interior models with the amount of 
C measured by the Cassini spacecraft requires that significant differential sedimentation of water and volatile 
species have taken place in Saturn's interior during its lifetime. 
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1. Introduction 

The standard giant planet formation scenario, namely 
the so-called core accretion model (Pollack ct al. 1996), 
has been recently extended by Alibert et al. (2004; 
2005a, b) to include migration of the growing planet and 
proto-planetary disk evolution. These calculations show 
that the extended core accretion model can lead to the 
formation of two giant planets closely resembling Jupiter 
and Saturn in less than 3 Myr (Alibert et al. 2005b, here- 
after A05b). Based on these formation models, Alibert 
et al. (2005c, hereafter A05c) also demonstrated that the 
volatiles enrichments measured in Jupiter's atmosphere 
can be explained in a way consistent with the constraints 
set from internal structure modeling (Saumon & Guillot 
2004, hereafter SG04). In this work, we follow the strategy 
adopted by A05c and utilize the thermodynamical condi- 
tions of the protoplanetary disk that led to the formation 
of Saturn (A05b) to calculate the volatile enrichments in 
its atmosphere. We then compare the calculated amount 
of required ices in Saturn with the content in heavy ele- 
ments predicted by the internal structure models of SG04. 



lations, we adopt two extreme values of C abundance that 
cover the domain of measurements published in the last 
25 years (see Tabled . From the use of a haze model to de- 
rive volume mixing ratios in the near-IR bands, C has been 
shown to be enriched by a factor of 3.4 ±0.9 (Kerola et al. 
1997; Kuiper Airborne Observatory (KAO)) in Saturn's 
atmosphere compared to its solar value 1 (Lodders 2003; 
see Table On the other hand, recent Cassini CIRS far- 
IR measurements lead Flasar et al. (2005) to infer a C 
enrichment of 8.8 ± 1.7 compared to solar from the fit of 
CH-4 rotational lines. 

Moreover, from an indirect determination, nitrogen 
has been shown to be enriched by a factor of 3.1 ± 0.6 
(Briggs & Sackett 1989; ground-based radio wavelength 
observations) in Saturn's atmosphere compared to solar. 
Briggs & Sacketts (1989) also estimated that S is enriched 
by a factor of ~ 12 compared to solar in Saturn. However, 
this latter estimation seems to be uncertain due to the 
difficulty of distinguishing in the microwave spectrum be- 
tween the contribution to the opacity of H2S from that of 
other compounds (Hersant et al. 2004). For this reason, 



The carbon abundance in Saturn's atmosphere has 
been the subject of a large number of studies. In our calcu- 
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1 Note that we use the central values of species abundances 
derived from Table |2] of Lodders (2003) in our calculations, 
contrary to A05b and A05c who adopted a set of values within 
the error bars to obtain the best fit possible of the volatile 
enrichments in Jupiter. 
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Fig. 1. Stability curves of the condensates considered in 
the present work, and evolutionary tracks of the nebula at 
12 and 9.8 AU. Abundances of various elements are solar, 
with C0 2 :CO:CH 4 = 30:10:1 and N 2 :NH 3 = 1 in vapor 
phase. The condensation curve of CO2 pure condensate 
(solid line) is shown together with that of the correspond- 
ing clathratc hydrate (dashed line). Species remain in the 
vapor phase as long as they stay in the domains located 
above the curves of stability. 

we assume in the present work that the abundance of S is 
still unknown. 

2. Formation of Saturn 

The model of Saturn's formation considered here is the 
one presented by A05b. In this model, Saturn forms from 
an embryo originally located at ~ 12 AU. The proto- 
planet migrates inwards by following proto-Jupiter's trail 
and stops at the current position of Saturn when the disk 
disappears. Planetesimals are accreted during the whole 
formation process, along proto-Saturn's migration path- 
way (between 12 and 9.8 AU). The resulting planet ex- 
hibits an internal structure comparable to that of the ac- 
tual Saturn described by SG04 2 , with a final core mass 
of ~ 6 M 8 and a total mass of accreted planetesimals of 
~ 13.2 M®. Finally, the thermodynamical conditions in 
the disk at the actual position of Saturn are used do calcu- 
late the composition of planetesimals that were accreted 
by the planet. 

3. Composition of planetesimals in Saturn's 
feeding zone 

In this work, we consider that abundances of all elements 
are solar and that O, C, and N exist only under the form 

2 Detailed internal structure models of Saturn developed by 
SG04 indicate a maximum mass of heavy elements present in 
the envelope between nearly and 10 M®, whereas the mass 
of the core varies between 8 and 25 M®. However, note that 
the minimum mass of the solid core might be decreased by up 
to ~ 7 depending upon the extent of sedimented helium. 



of H 2 0, C0 2 , CO, CH 4 , N 2 , and NH 3 in the solar nebula 
gas-phase. Concerning the initial ratios of C02:CO:CH4 
and N2:NH3 in the solar nebula, for the sake of consis- 
tency, we adopt the same ranges of values as the ones de- 
rived by A05c to fit the observed enrichments in volatiles 
in Jupiter. Hence, we assume C02:CO:CH4 = 40:10:1 to 
10:10:1 for N 2 :NH 3 = 1 and C0 2 :CO:CH 4 = 10:10:1 for 
N2:NH3 = 10 in the vapor phase of the solar nebula. The 
considered values of C02:CO:CH4 are consistent with ISM 
measurements (Allamandola et al. 1999; Gibb et al. 2004). 
The adopted values of N2:NH3 are compatible with ther- 
mochemical calculations in the solar nebula that take into 
account catalytic effects of Fe grains on the kinetics of N2 
to NH 3 conversion (Fegley 2000). In addition, S is assumed 
to exist under the form of H 2 S (H 2 S/H 2 = 0.7 x (S/H 2 ) Q , 
following A05c) and other sulfur compounds in the solar 
nebula (Pasek et al. 2005). 

Finally, we also assume that volatiles have been 
trapped during the cooling of the nebula in planetesimals 
either in form of pure condensates, as it is the case for 
CO2, or in form of hydrates or clathrate hydrates, as it is 
the case for NH 3 , N 2 , CH 4 , CO, Xe, Ar and Kr (see A05c 
for details). Note that the amount of available water ice is 
assumed to be sufficient to trap all volatiles except CO2. 
Once condensed, ices are assumed to decouple from the gas 
and to be incorporated into growing planetesimals which 
may subsequently be accreted by the forming Saturn. 

Figure n represents the cooling curves of the nebula at 
9.8 and 12 AU derived from the solar nebula model used 
by A05b, as well as the condensation curves for the var- 
ious ices considered in this work. The stability curves of 
clathrate hydrates are derived from Lunine & Stevenson 
(1985) whereas we used a fit to the experimental data 
for the pure CO2 condensate (Lide 1999). Figure ^ pro- 
vides the condensation sequence of the different volatiles 
initially existing in vapor phase inside Saturn's feeding 
zone. The intersection between the cooling curve and the 
stability curve of the different condensates also gives the 
thermodynamical conditions at which the different ices are 
formed. 

Using Eqs. 1 and 2 given in A05c, we can calculate 
1) the minimum abundance of water required to trap all 
volatiles (except CO2) in the solar nebula gas-phase to 
form clathrate hydrates or hydrates and 2) the composi- 
tion of ices incorporated in planetesimals that were ac- 
creted by proto-Saturn. Since, as shown in Table |31 the 
composition of planetesimals does not vary significantly 
along the migration path of Saturn, we assume that all 
the ices contained in accreted planetesimals have an iden- 
tical composition. 

4. Enrichments in volatiles in Saturn 

From the determination of the composition of the ices in- 
corporated in the planetesimals which depend on the ini- 
tial gas-phase composition in the solar nebula and the 
adjustment of the amount of heavy elements accreted by 
proto-Saturn, it is possible to calculate the expected en- 
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Table 2. Gas phase abundances (molar mixing ratio 
with respect to H2) of major species in the solar neb- 
ula (from Lodders 2003) for C0 2 :CO:CH 4 = 30:10:1 and 
N 2 :NH 3 = 1. 



Species i 


Xi 




Species i 


Xi 







1.16 x 10" 


-3 


N 2 


5.33 x 10" 


-5 


C 


5.82 x 10" 


-4 


NH 3 


5.33 x 10" 


-5 


N 


1.60 x 10" 


-4 


S 


3.66 x 10" 


-5 


H 2 


1.66 x 10" 


-4 


Ar 


8.43 x 10" 


-6 


C0 2 


4.26 x 10" 


-4 


Kr 


4.54 x 10" 


-9 


CO 


1.42 x 10" 


-4 


Xe 


4.44 x 10" 


10 


CH4 


1.42 x 10" 


-5 









Table 3. Calculations of the ratios of trapped masses of 
volatiles to the mass of H2O ice in planetesimals formed at 
9.8 and 12 AU in the solar nebula. Gas-phase abundance 
of H 2 is equal to 8.11 x 10" 4 (~ 4.9 times the value 
quoted in Table see text), and gas-phase abundances 
of elements are assumed to be solar (Lodders 2003) with 
C0 2 :CO:CH 4 = 30:10:1 and with N 2 :NH 3 = 1 in vapor 
phase in the solar nebula. The abundance of H 2 S is sub- 
solar (see text). 



9.8 AU 12 AU 



C0 2 :H 2 


8.91 x 10" 


-1 


9.13 x 10" 


-1 


CO:H 2 


1.46 x 10" 


-1 


1.45 x 10" 


-1 


CH 4 :H 2 


9.14 x 10" 


-3 


9.23 x 10" 


■3 


N 2 :H 2 


5.26 x 10" 


-2 


5.34 x 10" 


-2 


NH 3 :H 2 


4.67 x 10" 


-2 


4.66 x 10" 


-2 


H 2 S:H 2 


4.28 x 10" 


-2 


4.32 x 10" 


-2 


Ar:H 2 


9.74 x 10" 


-3 


9.94 x 10" 


■3 


Kr:H 2 


1.34 x 10" 


-5 


1.34 x 10" 


-5 


Xe:H 2 


2.44 x 10" 


-6 


2.45 x 10" 


-6 



richments in volatiles in the giant planet's atmosphere. 
In Table 0] we present the results of these computa- 
tions and compare them to the observed enrichments for 
C0 2 :CO:CH 4 = 30:10:1 and N 2 :NH 3 = 1 (our nominal 
vapor phase conditions) in the solar nebula gas-phase. 
It can be seen that our model reproduces the C and N 
ground-based enrichment measurements. The minimum fit 
of the airborne C and ground-based N measurements cor- 
responds to an O enrichment of 5.4 times the solar value. 
From this, we predict that S, Ar, Kr and Xe are over- 
abundant compared to solar by a factor of 2.1, 2.0, 2.2 
and 2.6 respectively. In addition, these enrichments result 
from the accretion of 5.6 M s of ices including 2.5 M® of 
water in the envelope of Saturn. The resulting total con- 
tent in heavy elements in Saturn is therefore compatible 
with the formation model of A05b and consequently with 
the internal structure model of SG04, provided the ratio 
of ices to rocks (I/R) is greater or equal to ~ 0.7. Table 
summarizes the ranges of required masses of ices in the 
envelope of Saturn, and the resulting O enrichments, to 



3 

match both C and N ground-based measurements if the 
initial gas phase conditions are varied. It can be seen that 
O is enriched at least between 5.1 and 7.6 times the so- 
lar abundance in Saturn's atmosphere while the minimum 
mass of incorporated ices varies between 5.4 and 7.6 M® 
in the giant planet. Under these conditions, the minimum 
I/R ratio required in Saturn varies between ~ 0.7 and 1.4 
in order to be compatible both with the formation model 
of A05b and the internal structure model of SG04. 

On the other hand, fitting the recent C measurement 
derived from Cassini data does not allow to retrieve a 
value of N compatible with the measured one. In the nom- 
inal case, fitting the minimum value of the Cassini mea- 
surement implies overabundances compared to solar of N, 
S, Ar, Kr and Xe by factors of 6.5, 5.5, 5.1, 5.6 and 6.7 
respectively. These calculations translate to an O enrich- 
ment of at least 13.9 times the solar value in the planet's 
envelope and require the accretion of at least 14.4 M® of 
ices including 6.6 M® of water. Following Table [SI if the 
initial gas-phase conditions are modified, the minimum 
mass of accreted ices varies between 13.7 and 18.3 M®, 
whereas the corresponding O enrichment ranges between 
13.1 and 18.3 times the solar value in the giant planet's 
atmosphere. In all these cases, the minimum amount of 
required ices exceeds the mass of heavy elements accreted 
by the Saturn's formation model of A05b. Moreover, this 
minimum mass also exceeds the mass of heavy elements 
present in the planet's atmosphere derived by SG04 (of 
the order of ~ 10M®). However, it is still compatible with 
the global amount of heavy elements predicted by SG04 
(< 30 Mg,). 

Table 4. Observed C and N enrichments in Saturn, and 
calculated enrichments in volatiles in the nominal model 
(C0 2 :CO:CH 4 = 30:10:1 and N 2 :NH 3 = 1). The observed 
values are taken from Kerola et al. (1997) (a), Flasar et al. 
(2005) (b) and Briggs and Sackett (1989) (c). Calculated 
enrichments are calibrated on the observed C enrichment 
owning the same column label. 



Species 


Observed 


Calculated 




(1) (2) 


(1) (2) 


C 


3.4 ± 0.9 a 8.8 ±1.7 b 


2.8 7.1 


N 


3.1 ±0.6 c 


2.5 6.5 


S 




2.1 5.5 


Ar 




2.0 5.1 


Kr 




2.2 5.6 


Xe 




2.6 6.7 



5. Discussion and conclusions 

In the framework of the clathrate hydrates trapping the- 
ory and using the gas-phase abundances of major volatile 
species in the solar nebula used by A05c to fit the volatiles 
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enrichments in Jupiter, we have calculated the enrich- 
ments in O, C, N, S, Xe, Ar, and Kr with respect to their 
solar abundances in the atmosphere of Saturn. Two differ- 
ent measurements of the abundance of carbon in Saturn 
have been considered for the calibration of our calcula- 
tions, namely the ground-based determination derived by 
Kerola et al. (1997) and the recent Cassini measurement 
as published by Flasar et al. (2005). 

Using the near-IR determination of C by Kerola et al. 
(1997), we are able to fit the measured value of N enrich- 
ment. Moreover, the predicted amount of ices accreted by 
Saturn is consistent with the formation model of A05b 
and the internal structure model of SG04, provided the 
minimum I/R ratio required in the giant planet varies 
between ~ 0.7 and 1.4, a value compatible with the re- 
cent Cassini determination of Phoebe's low density (1630 
Kg.m~ 3 ; Porco et al. 2005). 

On the other hand, using the recent far-IR determina- 
tion of C reported by Flasar et al. (2005), the measured 
value of N enrichment could no longer be reproduced. In 
addition, the amount of ices accreted by Saturn exceeds 
the mass of heavy elements predicted by the A05b model. 
This difference is likely to be due to the assumptions made 
by A05b in the calculation of the planetesimaPs disk (see 
also Alibert et al. 2005a for details). 

Furthermore, the amount of heavy elements required 
to explain the Cassini measurement that should be 
currently in Saturn's envelope appears incompatible with 
the results of the internal structure modeling by SG04 
(Table 0. In addition, in order to explain only the value 
of C enrichment given by Flasar et al. (2005), one needs 
to invoke the accretion of at least 10.8 of ices (includ- 
ing a minimum of 3.9 M® of water) in the envelope of 
Saturn with our nominal gas-phase conditions, an amount 
slightly exceeding the maximum mass of heavy elements 
predicted in this zone by SG04. The discrepancy is higher 
if one adds the amount of ices required to trap N2 and 
NH3 as clathrate hydrates or hydrates in planetesimals: 
the measured abondance of N by Briggs & Sackett (1989) 
requires as a minimum the additional accretion of 0.9 
Mg of ices, whereas assuming an enrichment of N of 6.5 
(as predicted by our calculations) requires the accretion 
of a minimum of 2.3 M® of ices. Therefore, in order to 
reconcile the Saturn's internal structure models and the 
C measurement made by the Cassini spacecraft, one must 
then suppose a sedimentation of some of the accreted 
water onto the core, although the major part of other 
accreted volatiles would have remained in the atmosphere. 
This assumption is compatible with the work of Fortney 
and Hubbard (2003; 2004) who argued that a significant 
oxygen depletion may have occurred in Saturn's atmo- 
sphere during its thermal history. This may be due to 
the immiscibility of oxygen in hydrogen at high tempera- 
ture and pressure conditions (Fortney and Hubbard 2003). 

Finally, a number of ground-based C measurements 
such as those of Buriez & de Bergh (1981) and Trafton 
(1985) are in agreement with the Cassini CIRS data, 



along with the Voyager IRIS data utilized by Courtin 
et al. (1984) and the ISO-SWS determination by Lellouch 
et al. (2001) (see Table EJ. This implies that the Cassini 
determination would appear more representative of the 
carbon abundance in Saturn's atmosphere than that of 
Kerola et al. (1997). On this basis, we note that our global 
oxygen abundance in Saturn is at least ~ twice higher 
than that predicted by Visscher & Fegley (2005) from net 
thermochemical reaction calculations in its atmosphere. 
On the other hand, the calculations of Visscher & Fegley 
(2005) are consistent with the actual Saturn's atmospheric 
composition and may not be valid for the whole planet. 
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Table 1. List of measurements of CH4 abundances in Saturn's atmosphere published since 1981. 



CH4/H2 (x 10 3 ) 


Solar enrichment 


Reference 


Measurement 


2 ± 0.5 


3.4 ± 0.9 


Kerola et al. (1997) 


KAO 


2 +i.o 

— 0.8 


— 1.4 


Encrenaz & Combes (1982) 


ground-based 


2 2 +0,8 

* —0.2 


Q + 1-4 
— U . d 


Killen (1988) 


ground-based 


2.5 


4.3 


Tomasko & Doose (1984) 


Pioneer IPP 


3 ± 0.6 


5.2 ± 1.0 


Karkoschka & Tomasko (1992) 


ground-based 


4 ± 2 


6.9 ± 3.4 


Buriez & de Bergh (1981) 


ground-based 


4.2 ± 0.4 


7.2 ± 0.7 


Trafton (1985) 


ground-based 


4.4 ± 1.2 


7.6 ± 2.1 


Lellouch et al. (2001) 


ISO-SWS 




7 7 +4.i 
'■'-3.3 


Courtin et al. (1984) 


Voyager IRIS 


5.1 ± 1.0 


8.8 ± 1.7 


Flasar et al. (2005) 


Cassini CIRS 



Table 5. Minimum mass of accreted water (M wat0 r), minimum mass of accreted ices (Mi ces ) and corresponding 
minimum O/H abundance (compared to solar value) in Saturn's atmosphere required to fit the observed C enrichments. 
Calculations are presented for the minimum fits of both KAO C and ground-based N measurements (Kerola et al. 
1997; Briggs & Sackett 1989) and Cassini C measurement (Flasar et al. 2005) in Saturn's atmosphere. Ranges of 
C02:CO:CH4 and N 2 :NH 3 gas-phase ratios used here are those determined by A05c (see text). 









N 2 


:NH 3 = 1 




N 2 :NH 3 = 10 






C0 2 :CO = 4 


C0 2 :CO = 


3 C0 2 :CO = 2 


C0 2 :CO = 1 


C0 2 :CO = 1 


KAO 


Mwatcr/Mffi 


2.2 


2.5 


3.0 


4.0 


4.6 




M ices /M® 


5.4 


5.6 


5.9 


6.6 


7.6 




O/H 


5.1 


5.4 


5.8 


6.6 


7.6 


Cassini 


M watcr /M ffi 


5.7 


6.6 


8.0 


11.0 


11.2 




M ices /M 9 


13.7 


14.4 


15.7 


18.2 


18.3 




O/H 


13.1 


13.9 


15.2 


18.0 


18.3 



